Reshaping elastic nanotubes via self-assembly of surface-adhesive nanoparticles 
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Elastic sheets with macroscopic dimensions are easy to deform by bending and stretching. Yet 
shaping nanometric sheets by mechanical manipulation is hard. Here we show that nanoparticle 
self-assembly could be used to this end. We demonstrate by Monte Carlo simulation that spherical 
nanoparticles adhering to the outer surface of an elastic nanotube can self-assemble into linear struc- 
tures as a result of curvature-mediated interactions. We find that nanoparticles arrange into rings or 
helices on stretchable nanotubes, and as axial strings on nanotubes with high rigidity to stretching. 
These self- assembled structures are inextricably linked to a variety of deformed nanotube profiles, 
which can be controlled by tuning the concentration of nanoparticles, the nanoparticle-nanotube 
diameter ratio and the elastic properties of the nanotube. Our results open the possibility of de- 
signing nanoparticle-laden tubular nanostructures with tailored shapes, for potential applications in 
materials science and nanomedicine. 



Direct mechanical manipulation can make macroscopic 
sheets conform to specific shapes. However, it is difficult 
to use mechanical means to reshape sheets with sizes in 
the micrometer range and smaller. An alternative at sub- 
micrometric scales would be the use of adhesive nanopar- 
ticles that induce local deformation and may catalyze 
global shape changes. Indeed, nanoparticles have been 
shown to drive the folding of graphene [1 , of thin films of 
silicon [2 and carbon nanotubes [3^ , the budding of fluid 
membranes by protein aggregation [H [5] , and the defor- 
mation of vesicles via the adhesion of nanoparticles |6| |7j . 

When nanoparticles adhere and deform a surface, effec- 
tive, curvature- mediated nanoparticle interactions arise 
as a result of the tendency of the surface to minimize the 
deformation caused by the nanoparticle imprints. One of 
the first studies of curvature- mediated interactions was 
that of Goulian and colleagues [8 , who calculated such 
effects in the context of protein aggregation in biolog- 
ical membranes. The effective pair interaction in fluid 
membranes is usually isotropic, has a Casimir-like func- 
tional dependence on particle separation [9^ and a non- 
trivial dependence on the protein shape. However, for 
elastic (tethered) surfaces, the overall effect of the forces 
at play is more complicated. Unlike fluid surfaces, which 
cannot withstand shear, elastic thin sheets can stretch 
in response to the forces applied by strongly adhering 
nanoparticles. If the nanoparticles are able to diffuse on 
the sheet, they should be able to self- assemble in a config- 
uration that reduces the mechanical cost of deforming the 
surface. However, the stretching energy associated with 
tethered surfaces imposes global geometric constraints to 
nanoparticle arrangements, and this leads to nontrivial 
many-body effects that extend across the surface. We 
expect the effective nanoparticle interactions to depend 
on the bending and stretching rigidities of the surface, 
its topology, and the relative location and specific extent 
of the local deformations imprinted by the nanoparticles. 
As a result, a significantly different — and phenomeno- 
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logically richer — behavior emerges when nanoparticles 
aggregate on elastic surfaces rather than in fluid inter- 
faces. Clearly, for multiple indenting (adhering) nanopar- 
ticles, one expects that the shape of the deformed sur- 
face is inextricably linked to the configuration of the ad- 
hered nanoparticles. Here, we use computer simulation 
to demonstrate that nanoparticles adhering and locally 
deforming the surface of an elastic nanotube can self- 
assemble into linear structures which themselves cause 
the global reshape of the nanotube. Our results suggest 
a new, nanoparticle-based route to design the effective 
elastic properties and the overall shape of a flexible sur- 
face. 

To this aim, we considered a simple, coarse-grained 
model of an elastic nanotube with bending and stretching 
deformation modes, and of adhesive spherical nanopar- 
ticles (smaller than the tube radius) which adhere to 
the nanotube via a generic short-range potential act- 
ing between any nanoparticle and the surface. Monte 
Carlo simulations of this model show that, for a wide 
range of nanoparticle-nanotube diameter ratios, nanopar- 
ticles arrange into one-dimensional strings. Moreover, 
the strings acquire a preferential orientation which de- 
pends on the elastic rigidity of the nanotube: on stretch- 
able nanotubes, strings of nanoparticles arrange as rings 
or helices; on nanotubes with high rigidity to stretch- 
ing they form axial strings. We analyzed these arrange- 
ments as a function of bending rigidity, nanoparticle area- 
density (that is, the relative area of the surface covered by 
nanoparticles) and the average nanoparticle bound area. 
We also find that the cross-sectional profile of the nan- 
otube can be shaped into ellipsoidal, triangular, rectan- 
gular and other regular and irregular forms, depending on 
the number and orientation of the self-assembled strings, 
the nanoparticle-nanotube diameter ratio and the elastic 
properties of the nanotube. Overall, our results suggest 
that nanoparticle organization and surface deformation 
can be controlled by tuning the concentration of adhe- 
sive nanoparticles and the mechanical properties of the 
surface. 

We modeled the elastic surface following a beads-and- 
string scheme \T0\. The nanotube consists of a triangu- 
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FIG. 1: Detail of the particle (left) and triangulated-surface 
(right) representations of the nanoparticle-nanotube system. 
Nanotubes are made of hexagonally arranged beads of diam- 
eter a. Nanoparticles have a diameter of lOcr and interact 
with the beads with energy per bead e when a bead is within 
a distance 2a from the nanoparticle surface. The nanoparti- 
cle bound area (bottom, shaded) is defined by the triangles 
whose all three vertices fall within the annular volume defined 
by the excluded volume of the particles and a sphere of radius 
6.5a. We report the nanoparticle bound area as the average 
value, Ab, for all nanoparticles bound to the nanotube and as 
a percentage of the nanoparticle surface. 



lated surface of hard beads connected by elastic links and 
arranged in a hexagonal lattice with an overall cylindri- 
cal shape, as sketched in Fig. [l] The triangulation of 
the network allows for an easy computation of elastic en- 
ergies, with rigidities to bending and stretching, hi and 
kf, respectively, characterizing the resistance to elastic 
deformation of the nanotube. The nanoparticles are de- 
scribed as hard spheres, and interact with the surface 
beads with a maximum energy per bead, e. Beads and 
nanoparticles do not interact with spheres of their own 
type, besides satisfying excluded volume constraints. We 
define the total energy of the model as the sum the Hel- 
frich bending elasticity, the stretching energy of links, 
and a nanoparticle-surface adhesive potential, which are 
written as 
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The discretized bending term has been used before [4|fTTl, 
and is the sum of the bending energies of all triangle pairs 
of the surface, r^j is the length of the link connecting 
vertices i and j; n^j and uji are the unit normal vectors 
to the triangles sharing that link; and Ej ^ij 
of the areas of the triangles sharing vertex i. We de- 
scribe the stretching contribution with a Finitely Exten- 
sible Nonlinear Elastic (FENE) potential [12], where Iq is 
the zero-energy length of a link and Im = 3<j its maximum 
length. We have chosen Iq = 1.229cr, which corresponds 
to an area fraction for the beads of 60% [13] • This means 
that a nanoparticle of lOa in diameter with = 5% is 



bound to 15 beads approximately. The adhesion term is 
a ramp well potential acting between each nanoparticle- 
bead pair, with a tangent distance At = 5.5<j. For conve- 
nience we chose A^ = 6.5cr and Xm = T.bcr. We have seen 
that the specific shape of the adhesive term of the po- 
tential does not affect qualitatively our results, provided 
that the potential is short-ranged [T4] . 

Clearly, both elastic constants and the adhesion en- 
ergy per bead affect the extent of the deformation. We 
find that a convenient geometric parameter to character- 
ize deformation is the average nanoparticle area that is 
bound to the nanotube (see caption of Fig. [l] for a precise 
definition). We simulate the nanotube with a periodic 
boundary in the direction parallel to its longitudinal axis 
X, and use periodic boundaries in the three cartesian co- 
ordinates for the nanoparticles. The initial configuration 
consists of an unperturbed nanotube, and of nanopar- 
ticles randomly placed in the simulation box. Monte 
Carlo moves at constant number of particles and tem- 
perature sample the configurational space by attempt- 
ing changes in the position of vertices and nanoparti- 
cles. In order to allow for stretch-free configurations, 
changes in the length of the box parallel to the nan- 
otube axis were also attempted, with constant pressure 
Px = 0. The results presented in this work correspond to 
simulations with nanoparticles of lOa in diameter. Ad- 
ditionally, we have simulated systems with larger, 15a- 
nanoparticles at various points in the parameter space, 
and have also compared our results to a similar simula- 
tion model [15]. System-size and surface-discretization 
effects do not appear to be significant. It is important to 
point out that the nanoparticle area-density affects self- 
assembly and deformation. Indeed, we observe that a 
nanoparticle area-density larger than ^ 0.2 is necessary 
for the nanoparticles to aggregate into linear structures. 
As the area-density increases, the linear structures be- 
come progressively interconnected, approaching the limit 
of homogeneous nanoparticle coverage. In this study we 
have focused in the intermediate area-density regime, for 
which considerable deformation and a larger degree of 
nanotube shapes can be accessed. 

The average nanoparticle bound area controls the de- 
gree of deformation of the nanotube. For 20 < < 
SO/cbT and < 5%, the elastic energy associated with 
the nanoparticle imprints is of the order of the thermal 
energy of the beads, and the nanoparticles can explore 
the entire surface on the nanotube. As increases, 
nanoparticles bound to the surface become less mobile. 
When 4% ^ < 10%, we observed that nanoparti- 
cles spontaneously organize into linear structures. For 
^ 10%, the adhesive energy needed to reach such 
bound area becomes larger than 10 /cbT per nanoparticle, 
and this makes the nanoparticle mobility on the surface 
too slow for the simulations to reach equilibrium within 
our simulation timeframe. This is schematically depicted 
in Fig. [2| where we also describe the self-assembled lin- 
ear structures as a function of stretching rigidity. For 
fairly stretchable nanotubes, nanoparticles self-assemble 
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FIG. 2: Nanoparticles self- assemble into ordered linear struc- 
tures at moderate adhesive strength, switching orientation 
from rings and helices to axial strings when the stretching 
rigidity is such that the stretching energy is approximately 
1 /cbT per bead. The boundaries of the simulated nanotubes 
match up to avoid boundary effects. 



into rings and helices, which can coexist on the same 
nanotube. For slightly less stretchable tubes, for which 
the stretching energy per bead is below ^ 1 /cbT, linear 
aggregates are shorter and their orientation can span all 
possible angles, although we have seen a preference for 
radial and axial orientations. Above that energy, we see 
preferentially axial arrangements. 

Axial strings are essentially stretch- free configurations. 
Indeed, a nanotube completely rigid to stretching cannot 
accommodate nanoparticle rings or helices because these 
configurations necessarily involve double-curvature im- 
prints, which have a stretching penalty. But the rigid 
nanotube can deform around axial strings by curving 
only in the radial direction. This is why axial strings 
only appear above 1 /cbT in stretching energy per bead. 
Indeed, there is a moderately sharp transition from the 
bending-dominated regime — for which rings and he- 
lices are the bending-minimizing configurations — to 
the stretching-dominated regime, in which the minimiza- 
tion of stretching energy becomes dominant. We be- 
lieve that the fact that the transition between the two 
regimes takes place within a relatively small range of 
stretching energies has its origin in the high relative 
cost of stretching with respect to bending. Indeed, the 
stretching-to-bending energy ratio due to an indentation 
of depth h in a thin shell of thickness t can be written as 
Es/E}) ~ (h/t)'^ ^ Kp/i^h? [16]. This is also applicable to 
a nanotube with radius R as long diS h <^ R. Clearly, for 
sufficiently thin surfaces, the stretching energy cost over- 



FIG. 3: a,b, cross-sectional cuts of deformed nanotubes. 
Symbols locate the corresponding energies in c. Nanopar- 
ticles have not been drawn for clarity, c, adhesion-stretching 
energies of configurations with axial strings (squares). One 
point corresponding to a helix (circle) is shown for compar- 
ison. Energies are given in units of feT per bead. The 
straight line is a fit to simulation data of axial strings with 
an (undeformed) nanotube-nanoparticle diameter ratio in the 
interval 1.017 — 2.033, with 15 — 55 nanoparticles, At in the 
5 — 8% range, and bending and stretching rigidities within 
20 - 100 feT and 120 - 960 feT/cr^ 



whelms that of the bending, effectively imposing a global 
constraint on the possible deformations of the nanotube. 

Next we discuss the shapes of the deformed nanotubes 
mediated by the linear aggregation of the nanoparticles 
bound to it, and also the connection between nanoparti- 
cle self-assembly and the shape of the nanotube. Fig. [3^ 
shows two nanotube profiles: one with a helically shaped 
string of nanoparticles and another with four, almost- 
parallel axial nanoparticle strings. Ring and helical struc- 
tures induce a screw-like type of deformation on the nan- 
otube. Axial strings, however, induce a larger variety of 
shapes, which can also show constrictions along the nan- 
otube axis if the strings are not fully parallel. Fig. ^ 
shows four different shapes corresponding to two (ellip- 
soidal profile), three (triangular profile) and four (square 
and bow-tie profiles) parallel axial strings. The more 
bent the nanotube profile is, the larger the adhesive en- 
ergy per bead of the configuration, as can be seen in 
Fig. This figure also shows that stretching the sur- 
face above 1 /cbT per bead comes at a very steep cost in 
adhesion. This is largely due to the fact that the sur- 
face wrapping an axial string does not form an entirely 
straight channel, but follows the curvature of the indi- 
vidual nanoparticles forming the string. 
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FIG. 4: Snapshots of nanotube profiles with high rigidity 
to stretching, for which nanoparticles self-assemble into axial 
strings. The strings have been constrained to span the length 
of the tube and to contain the same number of nanoparticles 
each. The profiles thus have a uniform shape in the axial 
direction. Profiles in the leftmost column have been drawn 
twice as big for clarity. For all profiles the nanoparticle radius 
is lOcr, At = 40 /cbT, kf = 500 k^Tja^ and = 10%. 



Fig. [3] suggests that, for the stretching-dominated 



regime, the shape of the nanotube can be controlled 
by the area-density of nanoparticles (or analogously, the 
number of axial strings) and the nanotube-nanoparticle 
diameter ratio. Fig. |4] illustrates this by showing snap- 
shots of the profiles obtained varying these two geometric 
parameters at constant bound area and elastic rigidities. 
We see that the profiles are fairly symmetrical but for 
the higher number of axial strings, for which the shape 
can presumably get trapped in metastable configurations. 
Also, a sufficiently high number of axial strings can cause 
a nanotube to collapse. 

By altering the profile of the nanotube, one has 
direct access to the effective bending rigidity of the 
nanoparticle-nanotube composite in the axial direction, 
as such rigidity depends linearly on the tube's cross- 
sectional moment of inertia. Furthermore, it should be 
possible to reversibly switch nanoparticle self-assembly 
on and off in experiments — for instance by tuning tem- 
perature, or by exploiting electrostatic or depletion inter- 
actions. Controlling nanotube shape may be relevant to 
applications of tubular nanomaterials [17h19), of bioac- 
tive nanotubes [20 , and in microfluidics [21] . We believe 
that experimental systems in which the here described 
coupling between nanoparticle self-assembly and nan- 
otube deformation occurs can be readily realized. Also, 
our approach may inspire alternative routes to manipu- 
lating the folding of thin films of silicon for photovoltaic 
applications [2 . 
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